LPT Orsay, 04-19 
RM3-TH/04-3 
ROMA- 1371/04 



The K —¥ 7T vector form factor 
at zero momentum transfer on the lattice a 



D. Becirevic 1 , G. Isidori 2 , V. Lubicz 3 ' 4 , G. Martinelli 5 , 
F. Mescia 2,3 , S. Simula 4 , C. Tarantino 3,4 , G. Villadoro 5 

1 Labor atoire de Physique Theorique, Universite Paris Sud, 
Centre d'Orsay, F-91405 Orsay-Cedex, France 

2 INFN, Laboratori Nazionali di Frascati, Via E. Fermi 40, 1-00044 Frascati, Italy 

3 Dip. di Fisica, Universita di Roma Tre, Via della Vasca Navale 84, 1-00146 Rome, Italy 

4 INFN, Sezione di Roma III, Via della Vasca Navale 84, 1-00146 Rome, Italy 

5 Dipartimento di Fisica, Universita di Roma "La Sapienza" , 
and INFN, Sezione di Roma, P.le A. Mow 2, 1-00185 Rome, Italy 



Abstract 

We present a quenched lattice study of the form factors f+(q 2 ) and fo(q 2 ) of the matrix 
elements (n\s'y f j i u\K). We focus on the second-order SU(3)-breaking quantity [1 — /+(0)], 
which is necessary to extract \V US \ from decays. For this quantity we show that 
it is possible to reach the percent precision which is the required one for a significant 
determination of | V us \ . The leading quenched chiral logarithms are corrected for by using 
analytic calculations in quenched chiral perturbation theory. Our final result, /r n (0) = 
0.960 ± 0.005 sta t ± 0.007 syst , where the systematic error does not include the residual 
quenched effects, is in good agreement with the estimate made by Leutwyler and Roos. 
A comparison with other non-lattice computations and the impact of our result on the 
extraction of \V US \ are also presented. 



a To appear in Nuclear Physics B. 



1 Introduction 



The most precise determination of the Cabibbo angle, or equivalently of the CKM matrix 
element \V US \ [1], is obtained from K — > -niv (K^) decays. The key observation which 
allows to reach a good theoretical control on these transitions is the Ademollo-Gatto theo- 
rem [2], which states that the form factors, f+(q 2 ) and fo(q 2 ), at zero four-momentum 
transfer, are renormalized only by terms of at least second order in the breaking of the 
SU(3) flavor symmetry. The estimate of these smallish corrections, i.e. of the difference 
of /+(0) = /o(0) from unity, is presently the dominant source of theoretical uncertainty 
in the extraction of \V US \. 

Chiral perturbation theory (CHPT) provides a natural and powerful tool to analyse 
the amount of SU(3) (and isospin) breaking due to light quark masses. As shown by 
Leutwyler and Roos [3], within CHPT one can perform a systematic expansion of the 
type /+(0) = 1 + h + h + ■ ■ where f n = 0[M K ^/(Anf^) n }. Because of the Ademollo- 
Gatto theorem, the first non-trivial term in the chiral expansion, f'2, does not receive 
contributions of local operators appearing in the effective theory and can be computed 
unambiguously in terms of M K , M w and f w (f 2 = —0.023 in the K° — > n~ case [3]). 
The problem of estimating /+(0) can thus be re-expressed as the problem of finding a 
prediction for 



This quantity is difficult to be evaluated since it depends on unknown coefficients of 0(p n ) 
chiral operators, with n > 6. Using a general parameterization of the SU(3) breaking 
structure of the pseudoscalar meson wave functions, Leutwyler and Roos estimated Af = 
(0.016 ± 0.008). Very recently, Bijnens and Talavera [4] showed that, in principle, the 
leading contribution to Af could be constrained by experimental data on the slope and 
curvature of fo{<l 2 )', however, the required level of experimental precision is far from the 
presently available one. For the time being we are therefore left with the Leutwyler- 
Roos result, and the large scale dependence of the 0(p 6 ) loop calculations [4, 5] seems to 
indicate that its 0.008 error might well be underestimated [6]. 

The theoretical error on \V US \ due to the Leutwyler- Roos estimate of /+(0) is already 
comparable with the present experimental uncertainty (see e.g. Ref. [7]). When the high- 
statistics K(2, results from KLOE and NA48 will be available, this theoretical error will 
become the dominant source of uncertainty on \V US \. Given this situation, it is then 
highly desirable to obtain independent estimates of /+(0) at the 1% level (or below). 
The purpose of the present work is to show that this precision can be achieved using 
lattice QCD. 

The strategy adopted in order to reach the challenging goal of a ~ 1% error, is based 
on the following three main steps: 

1. Evaluation of the scalar form factor fo(q 2 ) at q 2 = g^ ax = [Mk — M n ) 2 . 

Applying a method originally proposed in Ref. [8] to investigate heavy-light form 
factors, we extract /o(<?max) from the relation: 



Af = 1 + f 2 ~ /+(0) • 



(1) 



(<k\s 1q u\K)(K\u 1o s\h) = {Mk + M n ) 2 
(7r|M7 w|7r) (if Is70sl.fr) 4M K M n 



[f (q 2 max ;M K ,M n )} 2 , 



(2) 
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where all mesons are at rest. The double ratio and the kinematical configuration al- 
low to reduce most of the systematic uncertainties and to reach a statistical accuracy 
on /o(<7max) we U below 1%. 

2. Extrapolation of fo(q 2 nmy ) to /o(°)- 

By evaluating the slope of the scalar form factor, we extrapolate fo from g^ ax to 
q 2 = 0. We note that in order to obtain /o(0) at the percent level the precision 
required for the slope can be much lower, because it is possible to choose values of 
<?max ver Y close to g 2 = 0. 

For each set of quark masses we calculate two- and three-point correlation functions 
of mesons with various momenta in order to study the q 2 dependence of both fo(q 2 ) 
and f+{q 2 ). The latter turns out to be well determined on the lattice, whereas the 
former does not. We improve the precision in the extraction of fo(q 2 ) by constructing 
a new suitable double ratio which provides an accurate determination of the ratio 
fo{q 2 )/f+{q 2 )- We will define this ratio in Section 4. Fitting the g 2 -dependence of 
fo(q 2 ) with different functional forms, we finally extrapolate fo(q 2 nax ) to /o(0). The 
systematic error induced by this extrapolation, which is strongly reduced by the 
use of small values for g 2 nax , is estimated by the spread of the results obtained with 
different extrapolation functions. 

3. Subtraction of the leading chiral logs and chiral extrapolation. 

The Ademollo-Gatto theorem holds also within the quenched approximation [9], 
which has been adopted in this work. The leading 0(p A ) chiral corrections to /o(0), 
denoted by /| where the superscript q refers to the quenched approximation, are 
finite and can be computed unambiguously in terms of the 0{p 2 ) couplings of the 
quenched CHPT (qCHPT) Lagrangian [10, 11]. For these reasons, in order to get 
rid of some of the quenched artifacts we define the quantity 

R(M yr \ _ M _ l + f!(M K ,M n )-f (0;M K ,M^ 

x ' M " ] = (Alf " (AM^ ' (3) 

where AM 2 = M 2 K — M 2 and extrapolate it to the physical kaon and pion masses. 
The ratio R(Mx, M^): i) is finite in the SU(3)-symmetric limit; ii) does not depend 
on any subtraction scale; iii) is free from the dominant quenched chiral logs. We 
emphasize that the subtraction of f% in Eq. (3) does not imply necessarily a good 
convergence of (q)CHPT at order 0{p A ) for the meson masses used in our lattice 
simulations. The aim of this subtraction is to define the quantity Af in such a way 
that its chiral expansion starts at order 0(p 6 ) independently of the values of the 
meson masses. In the presence of sizable local contributions, we expect R(Mk, M n ) 
to have a smooth chiral behavior and to be closer to its unquenched analog than 
the SU(3)-breaking quantity [1 — /+(0)]. Extrapolating the values of R(M K , M n ) to 
the physical meson masses, we finally obtain 

Af = R(M^ S , M^ s ) x [(AM 2 ) 2 ]* = (0.017 ± 0.005 stat ± 0.007 syst ) , (4) 

where the systematic error does not include an estimate of quenched effects beyond 
0(p 4 ). Our result (4) is in good agreement with the estimate Af = (0.016 ± 0.008) 
obtained by Leutwyler and Roos in Ref. [3]. 
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The systematic error quoted in Eq. (4) is mainly due to the uncertainties resulting 
from the functional dependence of the scalar form factor on both q 2 and the meson 
masses. This error can be further reduced by using larger lattice volumes (leading 
to smaller lattice momenta) as well as smaller meson masses. In our estimate of 
Af discretization effects start at 0(a 2 ) and are also proportional to (m s — m u ) 2 , 
as the physical SU(3)-breaking effects. In other words, our result is not affected by 
the whole discretization error on the three-point correlation function, but only by 
its smaller SU(3)-breaking part. Discretization errors on Af are estimated to be 
few percent of the physical term, i.e. well within the systematic uncertainty quoted 
in Eq. (4). For a more refined estimate of these effects, calculations at different 
values of the lattice spacing are required. Finally, we stress again that the effects of 
quenching onto the terms of 0{p A ) are not estimated and thus not included in our 
final systematic error. 

Using the unquenched result for / 2 = —0.023 (in the K° — > n~ case) [3], our final estimate 
for f+ n (0) is given by 

ff^ (0) = 1 + / 2 - Af = 0.960 ± 0.005 stat ± 0.007 syst = 0.960 ± 0.009 . (5) 

The plan of the paper is as follows. In Section 2 we introduce the notation and give some 
details about the lattice simulation. Section 3 is devoted to the extraction of /o(<?max) 
by means of the double ratio method, while in Section 4 we study the q 2 dependence of 
the form factors and extrapolate the scalar form factor to q 2 = 0. The calculation and 
the subtraction of the quenched chiral logarithms as well as the extrapolation of Af to 
the physical masses is discussed in Section 5. The final estimate of /+(0), its comparison 
with non-lattice results and the impact on \V US \ are discussed in Section 6. Finally our 
conclusions are given in Section 7. 

2 Notations and lattice details 

The K — > n form factors of the weak vector current = s^^u are defined by 

(7c\p')\V,\K\p)) = Q [f l + {q 2 ){p + p')» + r_{q 2 ){p-p')»] , q 2 = (p - p'f , (6) 

where Cj is a Clebsch-Gordan coefficient, equal to 1 (2 -1 / 2 ) for neutral (charged) kaons. 
As usual, we express f!i(q 2 ) in terms of the so-called scalar form factor, 

im = /;(g 2 ) + M /_ Ml ni) ■ m 

By construction /q(0) = /+(0) and the differences between K° — > tt~ and K + — > 7r° 
channels are only due to isospin-breaking effects. In the following we shall concentrate on 
the K° — > 7r~ case and work in the isospin-symmetric limit b , dropping also the superscript 
% on the form factors. 

b For a discussion of isospin breaking effects see the recent review in Ref. [7] . 
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From Eq. (6) the form factors can be expressed as linear combinations of hadronic 
matrix elements. The latter can be obtained on the lattice by calculating two- and three- 
point correlation functions 

C^(t x ,t y ,p,p f ) = ^(O^yjV^x) O^ K (0))e-^ + ^'<^ , (8) 

x,y 

C K ^(t,p) = ^(^ W (t, a ? )O^ ) (0))e-^, (9) 

X 

where is the renormalized lattice vector current and 0\ = J75M, = (J75S are the 
operators interpolating n and K mesons. Since we do not consider SU(2)-breaking effects, 
we always use degenerate u and d quarks. 

Using the completeness relation and taking t x and (t y — t x ) large enough, one gets 



C^(t x ,t y ,p,p>) -j-^ ^^(n(p')\%\K(p)) e- E ^~ E ^-^ , (10) 

(t y -i,)->oo K TV 

C K ^(t,p(p')) ^ ^P^e- E ^\ (11) 

where E K = JW K + |p1 2 , E w = v^Wf + W and ^Z^~ ) = (0\O KM (0)\K(n)}. Then it 
follows 

C«*(t x ,t y ,p,p>) (n(p')\V,\K(p))_ _ 



c K (t x ,p) c*(t y -t x ,p') ^r^Zoo vz^z: 

Consequently the hadronic matrix elements {tt(p')\V /1 \K(p)) can be obtained from the 
plateaux of the l.h.s. of Eq. (12), once Z K and Z w are separately extracted from the 
large-time behavior of the two-point correlators (11). 

The procedure described above is the standard one to calculate form factors on the 
lattice. In this way, however, it is very hard to reach the percent level precision required 
for the present calculation. In the next Section we describe a very efficient procedure to 
get both the scalar form factor at q 2 = g^ ax and the ratio fo{q 2 ) / f+{q 2 ) with quite small 
statistical fluctuations. 

Numerical data have been obtained in the quenched approximation on a 24 3 x 56 
lattice, by using the plaquette gluon action at (5 — 6.2. In order to remove leading 
discretization effects, the non-perturbatively C(a)-improved Wilson action and fermionic 
bilinear operators are employed [12]. In case of the weak vector current one has 

v» = z v (1 + bv am * + am A ^ u + Cydu - sa ^ u) ; (13) 



where Zy is the vector renormalization constant, by and cy are 0(a)-improvement coef- 
ficients and the subscript t refers to the light u (or d) quark. 

We performed a first run by generating 230 gauge field configurations and choos- 
ing quark masses corresponding to four values of the hopping parameters, namely k e 



5 



{0.13390,0.13440,0.13490,0.13520}. Using K and tt mesons with quark content (k s k e ) 
and (keki) respectively, twelve different K — > n correlators (Cjf 7r ) have been computed, 
using both k s < kg and k s > kg, corresponding to the cases in which the kaon(pion) is 
heavier than the pion(kaon). In addition, using the same combinations of quark masses, 
also the three-point 7r — > K correlations (C™ K ) have been calculated. Finally, twelve 
non-degenerate K — > K and four degenerate n — > 7r three-point functions have been 
evaluated. 

In order to increase the number of different combinations of quark masses, we generated 
a second run of 230 independent gauge configurations with quark propagators computed 
at three new values of the hopping parameter, namely k G {0.13385,0.13433,0.13501}. 
In this way six different combinations of meson masses have been added to the first run, 
obtaining a total of 18 independent sets of pseudoscalar meson masses. 

The simulated quark masses are approximately in the range (0.5 4- 2) x m s , where 
m s is the strange quark mass, and correspond to K and it meson masses in the interval 
~ 0.5 4- 1 GeV. Though the simulated meson masses are larger than the physical ones, 
the corresponding values of g^ ax = (M K — M n ) 2 are taken as close as possible to q 2 = 
(see Table 2 in the next Section). 

As for the critical hopping parameter, we find k c = 0.135820(2) using the axial Ward 
identity. The physical values of meson masses in lattice units have been obtained using 
the method of "lattice physical planes" [13], in which the ratios Mk/Mk* and M n /M p 
are fixed to their experimental values. In this way we find 

[aM K ] phys = 0.189(2) , [aM^] phys = 0.0536(7) . (14) 

To improve the statistics, two- and three-point correlation functions have been aver- 
aged with respect to spatial rotation, parity and charge conjugation transformations. For 
the same reason we have chosen t y — T/2 in the three-point correlators, which allows 
to average the latter between the left and right halves of the lattice. Finally three-point 
correlation functions have been computed for the 10 different combinations of momenta 
(p,p') listed in Table 1. 



ap 


(0,0,0) 


(1,0,0) 


(1,1,0) 


(1,1,1) 


(0,0,0) 


(±1,0,0) 


(1,1,0) 


(0,1,1) 


(0,1,0) 


ap 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(0,0,0) 


(1,0,0) 


(1,0,0) 


(1,0,0) 


(1,0,0) 


(1,0,0) 



Table 1: Momentum configurations in units of27r/L = ix/12 of three-point correlators 
with p (p 1 ) being the kaon(pion) spatial momentum. 

The statistical errors are evaluated using the jacknife procedure for each run. To 
combine the results of a given observable obtained in the two runs, we have used a 
Monte Carlo procedure to generate values of the observable which are normally distributed 
around their averages with widths given by the jacknife errors. Such a procedure is 
adopted throughout this paper, including the chiral extrapolation of /o(0) described in 
Section 5. 
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• k s = 0.13390, k, = 0.13440 
■ k s = 0.13390, k, = 0.13490 
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4 8 12 16 20 24 28 

t/a 



Figure 1: Time dependence of the double ratio (15), averaged between the two halves of 
the lattice, for two sets of quark masses, given in the legend. The plateaux are taken from 
t/a = 12 to t/a = 16. 



3 Calculation of fo(q. 



2 

max 



Following a procedure originally proposed in Ref. [8] to study the heavy-light form factors, 
the scalar form factor has been calculated very efficiently at q 2 = g^ ax = (M K — M n ) 2 
(i.e. p = p' = q = 0) from the double ratio of three-point correlation functions with both 
mesons at rest: 

0,0) (7^(4,^,0,0) 

Cq^ (t X , ty, 0, 0) C^ltx, ty, 0, 0) 

When the vector current and the two interpolating fields are separated far enough from 
each other, the contribution of the ground states dominates, yielding 

+\ (7i\sjou\K) (K\uj s\ti) 2 2 (Mk + M n ) 2 

K {t x ,ty) — — = , -TT^TZT; ^. „.|_\ = i/oWmaxJj 777 77 ■ (lo) 



Ho\t x , t y ) = ^ KR ; ; ^ ^ • (15) 



The quality of the plateau for the double ratio can be appreciated by looking at Fig. 1. 

There are several crucial advantages in the use of the double ratio (15). First, there 
is a large reduction of statistical uncertainties, because fluctuations in the numerator 
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and the denominator are highly correlated. Second, the matrix elements of the meson 
sources, \JZk and \/~Z^ ) appearing in Eq. (10), cancel between numerator and denomina- 
tor. Third, the double ratio is independent from the improved renormalization constant 
Z v [l + by(am s + ame)/2] as well as from the improvement coefficient cy [see Eq. (13)] c . 
Therefore the knowledge of Zy, by and cy is not necessary to compute R and this ra- 
tio is automatically improved at order a 2 . Finally, the double ratio is equal to unity in 
the SU(3)-symmetric limit at all orders in the lattice spacing a. Thus the deviation of 
R from unity depends on the physical SU(3) breaking effects on /o(<?max) as we ^ as on 
discretization errors, which are at least of order a 2 (m s — mi) 2 (see Section 4.3). A similar 
consideration applies also to the quenching error, because the double ratio Ro(t x ,t y ) is 
correctly normalized to unity in the SU(3)-symmetric limit also in the quenched approx- 
imation. 

Having chosen t y = T/2, the three-point correlation functions are symmetric between 
the two halves < t x < T/2 and T/2 < t x < T. Statistical fluctuations, however, are 
independent in the two halves of the lattice. The best precision is reached when the 
double ratio is constructed in the two halves separately and then averaged. We obtain 
values of /o(<Zmax) with an uncertainty smaller than 0.1%, as it can be seen from Table 2 
and Fig. 2 (left). 



1.0125 



1.0075 



f (q m ax) 



1.0000 



r *. • 

by-*- 



-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 

a AM 



1.08 



1.04 



1.00 



0.96 



0.92 



0.88 



f + (q 2 ) 



-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 

a AM 



Figure 2: Values o// (9max) fl e ft) an d f+(q 2 ) (right) as a function of the SU (3) -breaking 
parameter a 2 AM 2 = a 2 M^ — a 2 M 2 . The dependence on g^ ax and q 2 is not explicitly 
shown in the plots. 

By replacing in Eq. (15) the time component of the vector current with the spatial 
ones and setting ap = ap' = 7r/12 (1, 0, 0) (the minimum non-zero momentum allowed on 

c The latter feature holds only when q — and the time component of the vector current is used. 
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k a -kt 


a 2 AM 2 


2 2 


/o(<7max) 


0.13390 -0.13440 


+0.01727(12) 


0.000737(10) 


0.99991(6) 


0.13390 -0.13490 


+0.03316(26) 


0.00365(6) 


1.0018(5) 


0.13390 -0.13520 


+0.04225(34) 


0.00745(15) 


1.010(2) 


0.13440 -0.13390 


-0.01777(12) 


0.000665(8) 


0.99974(6) 


0.13440 -0.13490 


+0.01634(15) 


0.000988(21) 


1.0008(2) 


0.13440 -0.13520 


+0.02565(23) 


0.00312(8) 


1.0057(9) 


0.13490 -0.13390 


-0.03496(30) 


0.00279(5) 


0.9993(4) 


0.13490 -0.13440 


-0.01674(17) 


0.000832(18) 


1.0002(2) 


0.13490 -0.13520 


+0.00948(10) 


0.000497(15) 


1.0014(2) 


0.13520 -0.13390 


-0.04491(46) 


0.00484(11) 


1.000(1) 


0.13520 -0.13440 


-0.02646(31) 


0.00222(6) 


1.0015(7) 


0.13520 -0.13490 


-0.00955(12) 


0.000416(13) 


1.0009(2) 


0.13385 -0.13433 


+0.01682(7) 


0.000669(6) 


0.99995(5) 


0.13385 -0.13501 


+0.03821(24) 


0.00519(6) 


1.0047(7) 


0.13433 -0.13385 


-0.01731(7) 


0.000610(6) 


0.99977(5) 


0.13433 -0.13501 


+0.02205(15) 


0.00193(3) 


1.0022(5) 


0.13501 -0.13385 


-0.04093(22) 


0.00380(5) 


0.9995(6) 


0.13501 -0.13433 


-0.02295(15) 


0.00154(2) 


1.0008(3) 



Table 2: Values of the hopping parameters k s and kt, a 2 AM 2 , a 2 q 2 iax and /o(?max) obtained 
with the double ratio method [see Eq. (16)]. 

our lattice), also the form factor / + can be extracted from the double ratio 

Cl^^x, ty, P, P) Cl K (t x , ty, P, J)) (- 2 )] 2 

C« K (t x ,ty,P,p)Cr(t X ,ty,P,p) i h ' ^ ] 

where f = (^fM^TW ~ \ f W+W) 

< ^ax- The uncertainty obtained in this 

case, however, is about twenty times larger than the one obtained for fo(q 2 nax ), as shown 
in Fig. 2. 

4 Momentum dependence of the form factors and ex- 
trapolation to q 2 = 

In this Section we perform the extrapolation of the scalar form factor from g 2 ^ to q 2 = 0. 

2 



To this end we need to evaluate the slope of /o, which in turn means to study the q 
dependence of the scalar form factor. We stress that, in order to obtain /o(0) at the 
percent level, the precision required for the slope can be much lower, since the values of 
<?max use d in our lattice calculations are quite close to q 2 = (see Table 2). We find indeed 
that a ~ 30% precision is enough for our purposes. 



9 



The form factors fo{q 2 ) and f+{q 2 ) can be expressed as linear combinations of the 
matrix elements (7r|V^|if). Their g 2 -dependence is obtained by studying these matrix 
elements determined according to Eq. (12) for the set of momenta listed in Table 1, using 
both time and spatial components of the vector current. 

The renormalization constant Zy of the lattice vector current and the improvement 
coefficients by and cy are needed in the calculation. In order to minimize the statistical 
fluctuations we have calculated Z v and by on the same sets of gauge configurations and 
combinations of quark masses used in the study of the form factors. 



4.1 Evaluation of Zy and by 

The vector current renormalization constant, Zy, and the improvement coefficient, by, 
can be extracted from the following relation 

- / 7 ' _^ = Z v (l + b v am q ) + 0(a 2 ) , (18) 

2C W (^T/2,0,0) 

where am q = (l/k s — l/k c )/2 for K — > K transitions and am q = (1/kt — l/k c )/2 for 
7r — > 7i ones. The results obtained for the l.h.s. of Eq. (18) show a clear linear dependence 
on the quark mass am q . A linear fit yields 

Z v = 0.7876 ± 0.0002 , by = 1.393 ± 0.006 , (19) 

in excellent agreement with the findings of Refs. [14]- [16]. In what follows we use the 
values of Zy and by given in Eq. (19) and adopt for the improvement coefficient cy the 
non-perturbative value c v = —0.09 from Ref. [14] d . 



4.2 g 2 -dependence of fo and / + 

The form factors fo(q 2 ) and f+(q 2 ) are shown in Fig. 3 as a function of q 2 . Note that, 
although we considered only ten independent combinations of meson momenta (see Table 
1), we have computed both K — > tt and ix — > K amplitudes, obtaining in this way twenty 
different values of q 2 . The form factor f+{q 2 ) is rather well determined with a statistical 
error of ~ 5 4- 20%, whereas for the scalar form factor the uncertainties turn out to be 
about 5 times larger. 

Our results for f+(q 2 ) are very well described by a pole-dominance fit, viz. 

as illustrated by the solid line in Fig. 3. The values obtained for the slope A + agree well 
with the inverse of the i^*-meson mass square for each combination of the simulated quark 
masses. A simple linear extrapolation in terms of the quark masses to the physical values 



d We note that the scalar form factor /o(<? 2 ) being proportional to the matrix element of d^V^ is 
independent of cy [see Eq. (13)]. The choice of cy affects only the determination of the form factor 
f+(q 2 ) and its numerical impact turns out to be smaller than the statistical uncertainty. 
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f + (q 2 ) f„(q 2 ) 




-0.30 -0.25 -0.20 -0.15 jO.10 -0.05 0.00 0.05 -0.30 -0.25 -0.20 -0.15 jO.10 -0.05 0.00 0.05 

(aq) (aq) 



Figure 3: Values of the form factors f+{q 2 ) (left) and fo(q 2 ) (right) as a function of the 
four-momentum square {aq) 2 , for the quark mass combination k s = 0.13390 and kg = 
0.13440. Open markers represent the results obtained with the double ratio method; the 
errors on these points are smaller than the size of the markers. In the left picture the 
solid line is the pole- dominance fit (20). 

yields A + = 0.026 ± 0.002 in units of the physical value of M 2 + , which is consistent with 
the PDG value A + = 0.028 ± 0.002 [17] e as well as with the recent measurement from 
KTeV A + = 0.02502 ± 0.00037 [18], obtained using a pole parameterization. The polar 
fit (20) also provides the value of /+(0). Due to the uncertainties in the determination 
of f+(q 2 ), however, the values obtained for /+(0) have statistical errors well above 1%, 
which do not allow to investigate SU(3)-breaking effects. The only way to keep these 
errors below 1% is to use the high-precision results obtained for io(<?max)- This in turn 
requires a drastic improvement in the evaluation of the scalar form factor as a function of 
q 2 . To this end, we have considered an alternative procedure based on the introduction 
of suitable double ratios, namely 

7~> / i . \ _ Cf n (t x , ty,P,p') CQ K (t X , ty,P,p') 

l[ x, y) ~ c**(t x ,ty,jr,F) c« K (t xl t yj p,p>) ' [ZL) 

from which a determination of the ratio of the form factors fo{q 2 ) / f+{q 2 ) can be obtained. 
The advantages of the double ratios (21) are similar to those already pointed out for 
the double ratio (15), namely: i) a large reduction of statistical fluctuations; ii) the 
independence of the renormalization constant Zy and the improvement coefficient by, 
and iii) Ri — > 1 in the SU(3)-symmetric limit. We stress that the introduction of the 
matrix elements of degenerate mesons in Eq. (21) is crucial to largely reduce statistical 

e Since experiments cover a narrow region of values of q 2 , the PDG value is obtained assuming a linear 
g 2 -dependence of the data. 
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fluctuations, because it compensates the different fluctuations of the matrix elements of 
the spatial and time components of the weak current. 

By denoting with Ri(q 2 ) the plateaux of the double ratios Ri(t x ,t y ) at enough large 
time distances, the ratio fo(q 2 )/f+(q 2 ) is given by 

fo(q 2 ) = 1 Q 2 ~(Ek + E' K ) (p + p') t + (E K + E n )(p + p') t R t (g 2 ) 

f + (q 2 ) M 2 K -Ml (E K + E> K )(p-p>) t -(E K -E n )(p + p>) z R^) ' 1 J 



where E' K = \/M\ + |p'| 2 . The g 2 -behavior of fo(q 2 ) can be then investigated after mul- 
tiplying the ratio fo(q 2 )/ f+(q 2 ), obtained from Eq. (22), by the values of f+{q 2 ) obtained 
with the standard procedure. The statistical uncertainties on Ri(q 2 ) and fo(q 2 ) turn out 
to be ~ 1 4- 20% and ~ 5 4- 20%, respectively. The quality of the results for fo(q 2 ), 
obtained in this way, is illustrated in Fig. 4 for one of the eighteen combinations of quark 
masses. These results can be directly compared with those obtained by using the standard 
method and shown in Fig. 3 (right). 

In order to extrapolate the scalar form factor to zero-momentum transfer we have 
considered three different possibilities, namely a polar, a linear and a quadratic fit: 

U{q 2 ) = /o^°(0)/(l-A ^V), (23) 
h{q 2 ) = rf^W-a + ArV), (24) 
fo(q 2 ) = ^(Ol.fl + Ar'^co? 4 ). (25) 

These fits are shown in Fig. 4. 

The polar, linear and quadratic forms fit our results for all quark masses with compa- 
rable x 2 values. The three fits provide both /o(0) and the slope A with values consistent 
with each other within the statistical uncertainties. Although tiny, the differences found 
for /o(0) are appreciable with respect to the size of the physical SU(3)-breaking effects. 
Consequently, for each combination of quark masses we consider all three extrapolated 
values of /o(0) and treat the difference as a systematic uncertainty in the rest of our 
analysis. The results obtained for /o(0) and the slope A are collected in Table 3. 

Before closing this Section we mention that our results for the slope A , extrapolated 
to the physical meson masses using a linear dependence in the quark masses, give in units 
of the physical value of M 2 + : \f oL} = 0.0122(22), A^ n ° = 0.0089(11) and Aj^ = 
0.0115(26). Our "polar" value \^ oL ^ is consistent with the recent determination from 
KTeV A = 0.01414 ± 0.00095 [18], obtained using a pole parameterization. 



4.3 Discretization effects 

Lattice artifacts on /o(0) due to the finiteness of the lattice spacing start at 0(a 2 ) and are 
proportional to (m s — mi) 2 , like the physical SU(3)-breaking effects. Indeed, the determi- 
nation of /o(<?max) i s affected only by discretization errors of 0[a 2 {m s — mi) 2 ], because the 
double ratio (15) is (9(a)-improved and symmetric with respect to the exchange m s <-> mi 
in the weak vertex. In addition, since q 2 nax is proportional to (m s — mi) 2 , 0(a 2 ) effects 
in the extrapolation from g^ ax to q 2 = also vanish quadratically in {m s — mi). Being 
in our calculation a~ x ~ 2.7 GeV, we expect discretization errors to be sensibly smaller 
than the physical SU(3)-breaking effects. 
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Figure 4: The form factor fo{q 2 ) obtained from the double ratios (21) for q 2 < g^ ax 
(full dots) and from the double ratio (15) at q 2 = g^ ax (open dot), for the quark mass 
combination k s = 0.13390 and kg. = 0.13440. The dot-dashed, dashed and solid lines 
correspond to the polar, linear and quadratic fits given in Eqs. (23-25), respectively. The 
inset is an enlargement of the region around q 2 = 0. 

A qualitative estimate of discretization errors can be obtained by considering the 
effect of the terms proportional to the improvement coefficient by, which cancel out at 
first order in the double ratio (15), but contribute at second order. A simple calculation 
shows that this correction is given by [1 + by (am s + ami) /2] 2 / '[(1 + b v am s ) ■ {1 + byam^)] ~ 
1 — a 2 by(m s — mi) 2 /A. The size of such a 0(a 2 ) contribution is a few percent of the whole 
result. Calculations performed at different values of the lattice spacing, combined with 
the extrapolation to the continuum limit, will certainly allow a quantitative estimate of 
discretization errors and a further reduction of this source of uncertainty. 

In Fig. 5 we plot the values of /o(0), obtained from the fit in q 2 given by Eq. (25), versus 
{a 2 AM 2 ) 2 . The results agree well with the quadratic dependence on a 2 AM 2 , expected 
from both physical and lattice artifact contributions, as shown in Fig. 5 by the solid line 
representing a naive fit to the form 

/ (0) = l-A(a 2 AM 2 ) 2 , (26) 

where A is a mass-independent parameter. 
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k s -kt 


/o (po,) (0) 


\(pol) 




,{lin.) 
A 




^ (quad.) 


0.13390 - 0.13440 


0.9981(3) 


2.5(3) 


0.9985(2) 


1.9(2) 


0.9983(3) 


2.2(4) 


0.13390 - 0.13490 


0.9912(16) 


2.9(4) 


0.9941(10) 


2.1(2) 


0.9915(15) 


2.8(4) 


0.13390 - 0.13520 


0.9858(41) 


3.2(6) 


0.9928(22) 


2.2(3) 


0.9858(40) 


3.2(6) 


0.13440 - 0.13390 


0.9981(2) 


2.4(3) 


0.9985(1) 


1.9(2) 


0.9984(2) 


2.1(3) 


0.13440 - 0.13490 


0.9975(5) 


3.3(5) 


0.9983(3) 


2.5(3) 


0.9976(5) 


3.2(5) 


0.13440 - 0.13520 


0.9940(21) 


3.7(7) 


0.9972(13) 


2.7(4) 


0.9946(19) 


3.5(7) 


0.13490 - 0.13390 


0.9923(11) 


2.5(3) 


0.9941(7) 


1.9(2) 


0.9931(11) 


2.3(3) 


0.13490 - 0.13440 


0.9978(4) 


2.9(4) 


0.9984(3) 


2.2(2) 


0.9980(4) 


2.7(4) 


0.13490 - 0.13520 


0.9993(4) 


4.3(9) 


0.9998(3) 


3.3(5) 


0.9996(4) 


3.6(9) 


0.13520 - 0.13390 


0.9882(22) 


2.5(3) 


0.9912(16) 


1.8(2) 


0.9894(21) 


2.2(3) 


0.13520 - 0.13440 


0.9951(13) 


2.9(4) 


0.9967(10) 


2.2(3) 


0.9955(11) 


2.7(4) 


0.13520 - 0.13490 


0.9993(3) 


3.9(7) 


0.9997(3) 


2.9(4) 


0.9995(3) 


3.4(7) 


0.13385 - 0.13433 


0.9982(3) 


2.6(4) 


0.9986(2) 


2.0(2) 


0.9986(3) 


1.9(5) 


0.13385 - 0.13501 


0.9863(34) 


3.5(7) 


0.9913(17) 


2.6(3) 


0.9923(34) 


2.4(7) 


0.13433 - 0.13385 


0.9981(2) 


2.7(4) 


0.9985(1) 


2.0(2) 


0.9985(3) 


2.1(4) 


0.13433 - 0.13501 


0.9945(15) 


4.0(8) 


0.9966(9) 


2.9(4) 


0.9971(14) 


2.7(7) 


0.13501 - 0.13385 


0.9880(16) 


3.0(4) 


0.9913(9) 


2.2(2) 


0.9902(16) 


2.5(4) 


0.13501 - 0.13433 


0.9959(8) 


3.2(5) 


0.9972(5) 


2.4(3) 


0.9969(8) 


2.5(5) 



Table 3: The form factor at zero-momentum /o(0) and its slope Ao, obtained from the 
three fits (23-25), for all combinations of the hopping parameters. 

5 Extraction of Af 

In order to determine the physical value of /o(0), we need to extrapolate our results of 
Table 3 to the physical kaon and pion masses. As discussed in the introduction, the 
problem of the chiral extrapolation is substantially simplified if we remove the effect of 
the leading chiral logs. To this purpose we consider the quantity 

A/ = 1 + h ~ /o(0) , (27) 

where fi represents the leading non-local contribution determined by pseudoscalar meson 
loops within CHPT. This procedure is well defined thanks to the finiteness of the leading 
non-local contributions both in the quenched and in the unquenched version of the effective 
theory. Having subtracted the leading logarithmic corrections, we expect Af to receive 
large contributions from local operators in the effective theory and, as in the case of the q 2 
slopes discussed in Section 4, to be better suited for a smooth polynomial extrapolation 
in the meson masses. 

As pointed out in the Introduction, the subtraction of f 2 in Eq. (27) does not imply 
necessarily a good convergence of CHPT at order 0(p 4 ) for the meson masses used in 
our lattice simulations. The aim of this subtraction is to define a quantity whose chiral 
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Figure 5: Values of /o(0), obtained from the quadratic fit (25), versus (a 2 AM 2 ) 2 . The 
solid line is the result of the linear fit (26). 

expansion starts at order 0(p G ) independently of the values of the meson masses, or a 
quantity better suited for a smooth chiral extrapolation. It is also worth to emphasise 
that the numerical impact of /f turns out to be very small (around or below 10 % with 
respect to the whole value of A/) for the highest meson masses used in our simulations. 
It increases (without exceeding 30 % of A/) in the case of lighter masses, where the chiral 
expansion should have a better convergence. 

5.1 Chiral loops within full QCD 

In the isospin-symmetric limit, within full QCD, the expression of the leading chiral 
correction f 2 is [3] 

3 3 

/2 = t^HttK + t:H v k , (28) 



where 



H 



PQ 



647T 2 / 2 



2 " Xi 2 

„ 2M 2 M 2 /M 2 \ 



(29) 



Note that f% is completely specified in terms of pseudoscalar meson masses and decay 
constants (f n ~ 132 MeV); it is negative (f 2 ~ —0.023 for physical masses), as implied 
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by unitarity [3, 19]; it vanishes as {M 2 K - M 2 ) 2 /(f 2 Ml) in the SU(3) limit, following the 
combined constraints of chiral symmetry and the Ademollo-Gatto theorem. 



5.2 Chiral loops within quenched QCD 

The structure of chiral logarithms appearing in Eqs. (28)-(29), is valid only in the full 
theory. In the quenched theory, the leading (unphysical) logarithms are instead those 
entering the one-loop functional of qCHPT [9, 10, 11]. We calculated this correction and 
present the results in this Section. Normalizing the lowest-order qCHPT Lagrangian as 
in Ref. [9], with a quadratic term for the singlet field <3> = str($) chosen as 

C%, = ^D,%D^ - , (30) 



we find 
where 



f! = H« K + H^ )K , (31) 



H q 



PK 



Ml 
96tt 2 /^ 



M 2 [M 2 K + M 2 P ) - 2aM 2 K M 2 P ( M\ 



{M 2 K -Ml) 2 l ° g \Ml ' " 



(32) 



with M 2 S -^ = 2Ml — M 2 . As anticipated, the one-loop result in Eq. (31) is finite because 
of the Ademollo-Gatto theorem, which is still valid in the quenched approximation [9] and 
forbids the appearance of contributions from local operators in f%. A proof beyond the 
one-loop level that the Ademollo-Gatto theorem (and more generally the Sirlin's relation 
[20]) holds within qCHPT can easily be obtained by applying the functional formalism 
to the demonstration in Ref. [20]. The latter needs only flavor symmetries which hold on 
the lattice also in the quenched case. 

It is worth to emphasize that the nature of the SU(3) breaking corrections in the 
quenched theory is completely different from that of full QCD: only contributions coming 
from the mixing with the flavor-singlet state are present and one finds /| > 0, which is 
a signal of the non-unitarity of the theory. For typical values of the singlet parameters 
(M ~ 0.6 GeV and a ~ [21]) and for the physical values of pion and kaon masses, one 
finds /| ~ +0.022. For the quark masses used in the simulation (and at the same values 
of M and a), the non-local correction is substantially smaller, /| < 0.003. This estimate 
is not changed significantly if a value of M as large as 1 GeV and/or a 20% variation of 
fa are considered. We anticipate that the effect of the above-mentioned variations of M 
and/or f n on our final result for Af is within 20% and it has been included in the quoted 
final uncertainties [see later Eq. (37)]. The reason why in Fig. 5 all values of /o(0) are 
found to be smaller than unity, as expected in a unitary theory, may be that for these sets 
of masses the contribution of local operators dominates over the quenched chiral logs. 

The high level of accuracy required in the study of the K i3 form factors can only be 
achieved once all possible sources of systematic errors are shown to be kept well under 
control. In a lattice calculation this includes in particular a study of finite volume effects. 
In our case such effects are expected to be smaller than other systematic effects, as we have 
explicitly checked by performing the analytical calculation presented in the Appendix and 
based on the techniques discussed in Refs. [22, 23]. 
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5.3 Extrapolation to the physical masses 



The evaluation of /| allows to express the lattice results in terms of the subtracted quantity 
Af defined in Eq. (27). To study the dependence of this quantity on the meson masses 
it is convenient to divide Af by (AM 2 ) 2 : 



R(M K ,M„ 



A/ 1 + f q 2 {M K ,M^)- f^-M K ,M n ) 

(a 2 AM 2 ) 2 (a 2 AM 2 ) 2 



(33) 



We have used /| computed from Eq. (31) at the corresponding values of meson masses, 
setting Mo = 0.6 GeV and a = 0.05. The error induced by a variation of Mq and a in 
a reasonable range of values is found to be negligible compared to the statistical error. 
We find that at the simulated masses the effect of /| does not exceed 30% of the value of 
[1 — /o(0)]. Since /| has a non-trivial, non-analytic dependence from the quark masses, 
we subtract its contribution at each value of the quark masses. 
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Figure 6: R(Mk, M n ) as a function of [a 2 M"^ + a 2 M 2 ] for the cases of the linear, quadratic 
and polar q 2 - extrapolation in the fit of fo(q 2 )- The solid lines represent the results of the 
linear fit (34)- Triangles indicate the values of R(Mk, M n ) extrapolated to the physical 
kaon and pion masses. 

As shown in Fig. 6, the dependence of R(Mk, M n ) on the squared meson masses is 
well described by a simple linear fit: 

R{M K , M n ) = c u + c 12 [{aM K ) 2 + (aM,) 2 ] , (34) 
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whereas the dependence on AM 2 is found to be negligible. We find an excellent \ 2 f° r a U 
the three sets of values of /o(0) obtained from the linear, quadratic and polar extrapolation 
in q 2 , see Fig. 6. In order to check the stability of the results, we have also performed 
quadratic and logarithmic fits, viz. 

R{Mk, M n ) = C21 + c 22 [(aM K ) 2 + (aM w ) 2 ] + c 23 [(aM K ) 2 + (aM n ) 2 } 2 , (35) 
R(M K ,M W ) = C31 + c 32 log[(aM x ) 2 + (aM w ) 2 ] . (36) 

It is reassuring to find that, while increasing the number of parameters in the chiral ex- 
trapolation as done in the quadratic fit (35) leads to larger uncertainties, the shift in 
the central values remains smaller than the errors. The logarithmic fit has been consid- 
ered in order to investigate the presence of possible important non-analytical terms in 
R(Mk, M 7r ) f . In Fig. 7 we show that linear, quadratic and logarithmic functional forms 
provide equally good fits of the data with consistent results for the extrapolation to the 
physical point. 





Linear Fit 


Quadratic Fit 


Logarithmic Fit 


Linear f (q 2 ) 


0.009 ±0.001 


0.008 ±0.003 


0.011 ±0.002 


Quadratic fo(q 2 ) 


0.013 ±0.002 


0.014 ±0.005 


0.018 ±0.003 


Polar f (q 2 ) 


0.015 ±0.002 


0.018 ±0.005 


0.020 ±0.003 





Linear Fit 


Quadratic Fit 


Logarithmic Fit 


Linear f (q 2 ) 


0.012 ±0.001 


0.013 ±0.003 


0.015 ±0.002 


Quadratic fo(q 2 ) 


0.017 ±0.002 


0.021 ±0.004 


0.023 ±0.002 


Polar f (q 2 ) 


0.019 ±0.002 


0.025 ±0.005 


0.026 ±0.003 



Table 4: Results for Af obtained from linear, quadratic and logarithmic fits of the ratio 
R(Mk, M n ) assuming linear, quadratic or polar functional forms for the extrapolation of 
the scalar form factor to q 2 = 0. Upper and lower tables correspond to different choices 
of the time interval chosen for the fits of the two-point correlation functions (see text). 

The extrapolated values of R(Mk, M n ) are converted into predictions for Af by using 
the values of the physical kaon and pion masses in lattice units given in Section 2. By this 
procedure, i.e. expressing both R(M K , M n ) and the meson masses in lattice units, we get 
rid of the error associated with the determination of the lattice spacing in the quenched 
approximation. In Table 4 we collect the values of Af obtained for all the functional 
forms assumed both in fitting the g 2 -dependence of fo(q 2 ) and in extrapolating the ratio 
R(M K , M w ) to the physical meson masses. We present two sets of results. The first set, 
given in the upper table, is obtained by determining the matrix elements y/Zx and \/~Zl 
[see Eq. (11)] from a fit of the two-point correlation functions in the same time interval 
chosen for the three-point correlators, namely t/a e [11,17]. For the second set, given 

We have also considered other fitting procedures, like a linear fit applied to the subset of data 
corresponding to [(aM^) 2 + (aM^) 2 ] < 0.18. The extrapolated value of R(Mk,M v ) at the physical 
masses is well within the spread of values obtained using the linear, quadratic and logarithmic fits. 
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Figure 7: Comparison among linear (34), quadratic (35) and logarithmic (36) fits of the 
ratio R(Mk, M n ) as a function of [a 2 M^ + a 2 M 2 ]. Triangle, square and diamond repre- 
sent the values of R(Mk, M n ) extrapolated to the physical meson masses. For illustrative 
purposes we have chosen the case in which a quadratic q 2 -dependence of f (q 2 ) has been 
used to extrapolate the scalar form factor to q 2 = 0. 



in the lower table, the same time interval chosen to extract the meson masses, i.e. t/a G 
[12,26], has been considered. 

From the spread of the values in Table 4 we quote our final estimate 

Af = R(Mf ys , Mf^ s ) x [(a 2 AI 2 ) 2 ]* 8 = (0.017 ± 0.005 stat ± 0.007 syst ) , (37) 

where the systematic error is dominated by the extrapolation of fo{q 2 ) in q 2 and of 
R(Mk, M n ) in the meson masses. 

The result in Eq. (37) is the main result of this paper. We have checked that all other 
sources of systematic errors, as those originated from different choices of the plateau 
intervals, from the uncertainties in the values of the renormalization constant Zy and 
of the improvement coefficients by and cy, and from the estimated size of discretization 
effects, are negligible with respect to the systematic error quoted in Eq. (37). 
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6 Estimate of /+(0) and impact on \V US \ 

The value of Af given in Eq. (37) is in excellent agreement with the estimate of Af made 
by Leutwyler and Roos using a general parameterization of the SU(3) breaking structure 
of the pseudoscalar- meson wave functions [3]. It must be stressed, however, that our 
result is the first computation of this quantity using a non-perturbative method based 
on QCD. Combining our estimate of Af with the physical value of f 2 from Eq. (28), we 
finally obtain 

ff*~ (0) = 0.960 ± 0.005 stat ± 0.007 syst (38) 

to be compared with the value f+° n (0) = 0.961 ± 0.008 given in Ref. [7] and quoted by 
the PDG [17]. 

6.1 Comparison with other approaches 

We now compare our results in Eqs. (37-38) with the recent evaluations of the 0{p & ) 
contributions to /+(0) made in Refs. [4, 6, 25]. In the chiral expansion the whole 0(p & ) 
term can be separated into two terms, coming respectively from meson loops and local 
terms. Whereas the whole 0(p 6 ) term is scale- and scheme-independent, the two terms 
separately depend on the renormalization scale and scheme. Leutwyler- Roos assumed that 
their quark-model estimate represents the whole physical 0(p 6 ) term. In Refs. [4, 6, 25] 
this result was interpreted instead to represent the contribution due to the local terms 
only, at the renormalization scale of the p-meson mass. The full 0(p 6 ) contribution was 
then obtained by adding the renormalized loop amplitude computed in Ref. [4]. In this 
way the whole 0{p & ) correction is found to be very small, though within a large error. The 
compensation between local and non-local terms at 0(p 6 ) strongly depends on the choice 
of the renormalization scale. Indeed, varying the renormalization scale from /x^jg ps 1 GeV 
to /ij^g ps 0.5 GeV, the loop contribution to /+(0) changes from +0.4% to +3.5% [6]. 

As pointed out in Ref. [25], the use of dispersion relations could help to solve the 
scale ambiguity. By means of dispersion relations one obtains complementary constraints 
on the slope and curvature of /o(<? 2 ) that, in turn, can be used in conjunction with the 
chiral calculation of Ref . [4] to determine /+(0). Present data, however, are not accurate 
enough: the estimate of 0(p 6 ) local terms given in Ref. [25] can only be interpreted as 
a consistency check of the Leutwyler-Roos ansatz, rather than as a truly independent 
determination. Finally, it must be stressed that all the purely analytic approaches to the 
estimate of /+(0) suffer from an intrinsic uncertainty due to the value of fx/ fn [26]. 

Both these ambiguities (the scale uncertainties and the dependence on fx/fir) are not 
present in our approach. By construction, Af is a scale-independent quantity which in- 
cludes all contributions beyond 0(p A ). We determine this quantity only in terms of meson 
masses. The main sources of uncertainty in our approach come from the (^-dependence of 
the form factor, the chiral extrapolation and the quenched approximation, but we argue 
that such uncertainties are included in the systematic error quoted in Eq. (37). 
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6.2 Impact on \V US \ 

Given the good agreement between our estimate of f+° n (0) and the value of Refs. [7, 17], 
our result implies a negligible shift in the determination of \V US \ from Kg 3 with respect to 
those works. Assuming the same experimental inputs, the value computed with our form 
factor, applying the updated radiative and isospin-breaking corrections of Refs. [6, 24], 
is \V US \ = 0.2202 ± 0.0025 to be compared with \V US \ = 0.2196 ± 0.0026 of Refs. [7, 17]. 
As pointed out by several authors, this estimate shows a m 2a deviation from the CKM 
unitarity relation, which yields \V US \ = 0.2278 ± 0.0020, taking into account the recent 
determination |T4<f| = 0.9737 ± 0.0007 from Ref. [27]. According to our analysis and 
assuming that the quenching effects do not drastically change our findings, this deviation 
cannot be attributed to an underestimate of the amount of SU(3) breaking in f+° n (0). 

We stress, however, that the above estimate of \V US \ takes into account only the (old) 
published data. Recently, new results have been presented [28, 29, 30, 31]. For 
instance, using only the high-statistics K e3 results of Refs. [28] and [29], we find sub- 
stantially higher values, namely \V US \ = 0.2275 ± 0.0030 and \V US \ = 0.2255 ± 0.0025, 
respectively, which are in good agreement with CKM unitarity. The KTeV findings [29] 
are also confirmed by the preliminary results presented by KLOE [30] and NA48 [31] at 
the ICHEP '04 Conference. 

7 Conclusions 

We presented a quenched lattice study of the K — > ir vector form factor at zero-momentum 
transfer. Our calculation is the first one obtained by using a non-perturbative method 
based only on QCD, except for the quenched approximation. Our main goal is the deter- 
mination of the second-order SU(3)-breaking quantity [1 — / + (0)], which is necessary to 
extract \V US \ from decays. In order to reach the required level of precision we employed 
the double ratio method originally proposed in Ref. [8] for the study of heavy-light form 
factors. We found that this approach allows to calculate the scalar form factor fo(q 2 ) at 
q 2 = g^ ax with a statistical uncertainty well below 1%. 

A second crucial ingredient is the extrapolation of the scalar form factor to q 2 = 0. 
This was performed by fitting accurate results obtained using suitable double ratios of 
three-point correlation functions. The systematic error due to finite volume effects was 
evaluated and discretization errors were qualitatively estimated to be a few percent of the 
deviation of /o(0) from unity. Calculations performed at different values of the lattice 
spacing, combined with the extrapolation to the continuum limit, will certainly allow a 
more quantitative estimate of discretization errors and a further reduction of this source 
of uncertainty. Nevertheless it is reasonable to conclude that the uncertainties coming 
from the functional dependencies of the scalar form factor on both q 2 and the meson 
masses are the dominant contribution to the systematic error. 

The leading chiral artifacts of the quenched approximation, /|, were corrected for 
by means of an analytic calculation in quenched chiral perturbation theory. After this 
subtraction, the lattice results have been smoothly extrapolated to the physical meson 
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masses. Our final value for /. 
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(0) is 



/- 



?K°7r- 
+ 



(0) = 0.960 ± 0.005 stat ± 0.007 syst = 0.960 ± 0.009 , 



(39) 



where the systematic error does not include an estimate of quenched effects beyond 0(p 4 ). 

The impact of our result on the determination of \V US \ was also addressed. Using the 
(old) published K E3 data from Ref. [17], we obtain \V US \ = 0.2202 ± 0.0025, which still 
implies a ~ 2a deviation from the CKM unitarity relation. We stress that, according 
to our analysis, such a deviation should not be attributed to an underestimate of the 
amount of SU(3) breaking effects in f+ n (0). Using only the recent high-statistics K e3 
results of Refs. [28, 29] one finds substantially higher values, \V US \ = 0.2275 ± 0.0030 and 
\V US \ = 0.2255 ± 0.0025, which are in good agreement with CKM unitarity. 

In order to reach a precision better than 1% in the determination of the Cabibbo angle 
both theoretical and experimental improvements are needed. As for the experimental side, 
new high-statistics results on both charged and neutral K i3 modes are called for. As for the 
theoretical side, the next important steps are: i) to remove the quenched approximation; 
ii) to decrease the values of the simulated meson masses in order to gain a better control 
over the chiral extrapolation of lattice results, and iii) to use larger lattice volumes for 
decreasing lattice momenta, in order to improve the determination of the slope of the 
scalar form factor. 
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Appendix 

In this Appendix we present the results of the analytic calculation of finite volume cor- 
rections to f'2, both in the quenched and unquenched theory. 

An important advantage of the effective theory approach is that it allows to evaluate 
(and eventually to correct for) the lattice artifacts due to a finite volume. This is simply 
achieved by imposing periodic boundary conditions to the wave functions of the chiral 
fields. This approximation represents a good description of the finite-volume effects as 
long as the size L of the lattice satisfies the condition L ■ M n ^ K 3> 1 [22], which is well 
satisfied in the present simulation. 

Defining the finite- volume shifts as 



fa = /2I1, - h , 



(40) 



we find that the corrections to Eqs. (28) and (31) are given by 




-2^{L,M V ,M K )-2^(L,M V ,M K )\ , 



(41) 
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(5M|a-4M>-M 2 ) 



M 2 K -Ml 

K 7T 



($i(L,M w )-ej(AM (rt) )) 
+3 (M 2 - a M 2 ) £ f (L, AQ + 3 (M 2 - a M? aS) ) £ f (L, M (s5) ) 

+8 (M 2 5) a - Ml) £s(L, M^-, M (m - } ) + 12 (M 2 a - M 2 ) Mr-, M ff ) 

2 2 

+ 12 (ML) a - Ml) &(L, Mjf, M (sg) )l , (42) 



in the unquenched and quenched cases, respectively. The functions £ s are defined as 



UL,M 1 )=UL,M 1 ,M 1 ) 
UL,M U M 2 ) 



1 



1 



dTT .-l e -M?r 



1 .3 



r 



(4^)3/2 



r s 



drr s 



_m(+m£ 



sinh < 



T 



x 



^(L,M 1 ,M 2 ) = 



9 



T 



(4ttt) 3 / 2 



;6(i,Mi,M 2 ), 



(43) 



where 6(x) = Xl^L-oo e ~ :cn2 , and they represent the differences between finite volume 
sums and infinite volume integrals (see Ref. [23] for more details). While S L represents 
a small correction in the range of masses used in the simulation, the effect of 8 q L is not 
negligible with respect to the size of the quenched logs in Eq. (31). Nevertheless, given 
the smallness of the leading chiral corrections at the simulated meson masses, we find 
that the inclusion of finite volume corrections has a negligible effect on our final estimate 
of Af. 
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